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Abstract — In this article, we propose a design method of self¬ 
interference cancelers for wireless relay stations taking account 
of the hasehand signal suhspace. The problem is first formulated 
as a sampled-data H°° control problem with a generalized sam¬ 
pler and a generalized hold, which can be reduced to a discrete¬ 
time -induced norm minimization problem. Taking account 
of the implementation of the generalized sampler and hold, we 
adopt the filter-sampler structure for the generalized sampler, 
and the uspampler-filter-hold structure for the generalized hold. 
Under these implementation constraints, we reformulate the 
problem as a standard discrete-time H°° control problem by 
using the discrete-time lifting technique. A simulation result is 
shown to Illustrate the effectiveness of the proposed method. 

I. INTRODUCTION 

In wireless communications, relay stations have been used 
to relay radio signals between radio stations that cannot 
directly communicate with each other due to the signal 
attenuation. More recently, relay stations are used to achieve 
spatial diversity called cooperative diversity to cope with 
fading channels [1]. On the other hand, it is important to 
efficiently utilize the scarce bandwidth due to the limitation 
of frequency resources [2], while conventional relay stations 
commonly use different wireless resources, such as fre¬ 
quency, time and code, for their reception and transmission 
of the signals. For this reason, a single-frequency full- 
duplex relay station, in which signals with the same carrier 
frequency are received and transmitted simultaneously, is 
considered as one of key technologies in the fifth generation 
(5G) mobile communications systems [3]. In order to realize 
such full-duplex relay stations, self-interference caused by 
coupling waves is the key issue [4]. 

Fig. [T] illustrates self-interference by coupling waves. In 
this figure, radio signals with carrier frequency / [Hz] are 
transmitted from the base station (denoted by BS). One 
terminal (denoted by Tl) directly receives the signal from 
the base station, but the other terminal (denoted by T2) is so 
far from the base station that they cannot communicate di¬ 
rectly. Therefore, a relay station (denoted by RS) is attached 
between them to relay radio signals. Then, radio signals 
generated by BS with the same carrier frequency / [Hz] 
are transmitted from RS to T2, but also they are fed back 
to the receiving antenna of RS directly or through reflection 
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Fig. 1. Self-interference 


objects. As a result, self-interference is caused in the relay 
station, which may deteriorate the quality of communication 
and, even worse, may destabilize the closed-loop system. 

To solve this problem, the authors have proposed a 
design method of self-interference cancelers that stabilize 
the closed-loop system and also suppress the continuous¬ 
time effects of self-interference, based on sampled-data 
control [5]. This method takes transmitted signals as Al¬ 
tered Lf signals and minimizes the norm (or the L^- 
induced norm) of the error system between received and 
transmitted signals. However, in many cases, signals used 
in communications belong to a subspace of Lf characterized 
by basis functions. This subspace is much smaller than the 
space of Altered Lf signals and hence the method becomes 
conservative. 

In this paper, we assume that we know the signal sub¬ 
space spanned by known basis functions, and utilize the 
characteristic of the subspace for self-interference canceler 
design. The problem of the canceler design is formulated by 
sampled-data H°° control problem with generalized sam¬ 
pler/hold. We show that the problem is reducible to a stan¬ 
dard discrete-time H°° control problem, which can be easily 
solved by numerical computation. We also show simulation 
results to illustrate the effectiveness of the proposed method. 

The remainder of this paper is organized as follows. In 
Section |II] we derive a mathematical model of a relay station 
considered in this study. In Section [III] we formulate the 
design problem as a sampled-data H°° control problem, 
which is transformed to a standard discrete-time H°° control 

















problem. In Section |V] an example is shown to illustrate the 
effectiveness of the proposed method. In Section|Vll we offer 
concluding remarks. 


Notation 


Throughout this paper, we use the following notation. 
We denote by the Lebesgue space consisting of all 
square integrable functions on [0, oo) endowed with the inner 
product 

poo 

{x,y) = / y{tY x{t)dt, 

and the norm ||a;|| = {x, x). We also denote by the set 

of all square summable sequences on {0, 1 , 2 ,.. .} endowed 
with the inner product 

OO 

{x,y)i-^ ='^y[k]^x[x] 

k=0 


and the norm ||x ||^2 = yj{x, x)p. The symbol t denotes the 
argument of time, s the argument of Laplace transform, and 
z the argument of Z transform. These symbols are used to 
indicate whether a signal or a system is of continuous-time 
or discrete-time. The operator with nonnegative real 

number L denotes a continuous-time delay operator and the 
operator z~^ with nonnegative integer / denotes a discrete¬ 
time shift operator. 

A transfer function of a discrete-time system with state- 
space matrices A, B, C, D is denoted by 


■ A 

B ' 

C 

D 


= C{zl - A)-^B + D. 


II. SYSTEM MODEL 


In this section, we derive a mathematical model of a relay 
station with coupling waves. 

There are mainly two types of relaying protocols in 
wireless communications [6]. One protocol demodulates 
received signals into digital signals, which are then ampli¬ 
fied and transmitted. This is called an amplify-and-forward 
(AE) protocol. The other protocol is called a decode-and- 
forward (DE) protocol, in which the received signals are 
demodulated, decoded and detected before transmission. The 
most essential difference between the AE protocol and the 
DE protocol is whether to detect received signals or not. 
In the DE protocol, the detection may reduce the likeli¬ 
hood of instability. However, the AE protocol requires quite 
lower implementation complexity than the DE protocol [6]. 
Throughout this study we assume that relay stations operate 
with the AE protocol. 

Fig. |2] depicts a typical AE single-frequency full-duplex 
relay station implemented with a digital canceler. A radio 
wave from a base station (BS) is accepted at the receiving 
antenna through Channel 1. Then, the radio frequency (RE) 
signal is demodulated to a baseband signal by the RE 
demodulator, and converted to a digital signal by the analog- 
to-digital converter (ADC). The obtained digital signal is 
then processed and the effect of loop-back interference 
mentioned in the previous section is suppressed by the digital 


Channel 2 



Fig. 2. AF single-frequency full-duplex wireless relay system 



Fig. 3. Mathematical model of the wireless relay system 


signal processor (DSP), which is converted to an analog 
signal by the digital-to-analog converter (DAC). Einally, the 
analog signal is modulated to an RE signal, amplified by a 
power amplifier (PA), and transmitted to terminals (T) with 
which we communicate by the transmission antenna through 
Channel 2. Inevitably, this signal is also fed back to the 
receiving antenna of the relay station. This is called coupling 
wave and causes loop-back interference, which deteriorates 
the communication quality. Practically, radio waves may vary 
because of noise and fading on Channel 1 and Channel 2. In 
this study, to concentrate the effect of suppressing coupling 
waves, we assume that the characteristics of the channel 1 
and 2 in Pig.|2] are modeled as the identity. It is notable that 
the assumption is not restrictive when the characteristic of 
the channels are known (e.g. channels affected by frequency- 
selective fading). In this case, we can include the channel 
model into the system model. 

Pig. [3 shows a simplified block diagram of the relay 
station. We model PA in Pig. [3 as a static gain denoted 
by a. The RE modulator is denoted by Ai and the RE 
demodulator by V. In this study, we adopt a multipath delay 
system, ^ channel model, where M is 

the multipath number, Vj > 0 is the gain, and Lj > 0 is 
the delay time of the j-th path. The block named “Digital 
System” includes ADC, DSP and DAC in Pig. |3 

In this article, we consider the quadrature amplitude 
modulation (QAM), which is used widely in digital com¬ 
munication systems, as a baseband modulation method. RE 
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Fig. 4. Digital system 


modulation of QAM transforms a transmission symbol into 
two orthogonal carrier waves, that is, a sine wave and a 
cosine wave. We assume that the transmission signal u{t), 
which is the output of “Digital System” in Fig. [3 is given 
by 


u{t) = ^ Ud[n]cj){t - nT), Ud[n] = 

n—0 


ui[n] 
U 2 [n] 


G 


where G is a basis function (also known as a pulse¬ 
shaping function), and T is the symbol period. The elements 
Ui[n], U2[n] are respectively the in-phase and the quadrature- 
phase components of the n-th transmission symbol. The 
same assumption is also applied to the received baseband 
signal from BS 

w = Vv. 


In this study, we define 

(t>n{t) = (t>{t - nT), n = 0,1,2,... 


and the set {^o, <^i, ■ ■ ■} is an orthonormal system in that 


is. 




1, 

0, otherwise. 


Note that this assumption is satisfied for a square root 
raised cosine (SRRC) pulse, which is widely used in digital 
communications [7]. 

Next we model the “Digital System” block in Fig. [3 as 
a sampler-filter-hold system shown in Fig. 0] In this block 
diagram, 5^ is a generalized sampler defined by 


50 : 


yid)] 1 

.2/2 (i)J / tg[o,oo) 


_{y2, <Pn)_ f n=0,l,2,... 


K{z) is a digital filter that we design for self-interference 
cancelation. is a generalized hold defined by 


: 


ui[n] 

U2[n] 




n=0,l,2,... 



ui[n] 

U2[n] 




t€[0,oo) 


Finally, after some transformations [5], we obtain a relay 
station model depicted in Fig.|3 which is equivalent to Fig. [3 
In this block diagram, we define 


M 

Qis) = 

i=i 

The input signal w is defined hy w = Vv as mentioned 
above, and we assume that w can be represented by a linear 
combination of .... By the block diagram, we can 

see that the relay station with self-interference is a feedback 


cos 2Tr fLj sin 27r/Lj 
— sin 271 fLj cos 27r fLj 



Fig. 5. Transformed wireless relay system 



Fig. 6. Sampled-data eri'or system 


system. In practice, the gain of the power amplifier (PA) 
in Fig. 13 is very high (e.g., a = 1000) and the loop gain 
becomes much larger than 1, and hence we should discuss the 
stability as well as self-interference cancelation. To achieve 
these requirements, we formulate the design problem as a 
sample-data H°° control problem. 

III. PROBLEM FORMULATION 

In this section, we formulate the design problem as a 
sampled-data H°° control problem based on the model 
discussed in the previous section. Then we show that the 
problem is reduced to a problem of discrete-time ^^-induced 
norm minimization. 

A. Design Problem via Sampled-Data H°° control 

The objective is to reduce the error between the input w 
and the output u shown in Fig. |3 In many cases, the output 
u is allowed to be delayed against the input w, and hence 
we consider the delayed error 

z(t) = w{t — mT) — u{t), 

where m is a non-negative integer. The corresponding error 
system is shown in Fig. |3 

Then, as mentioned above, the received signal w is rep¬ 
resented by a linear combination of fn- Let us define the 
signal subspace W by 

W = (1) 

where span{(/)o, fi ,..is the space consisting of all finite 
linear combinations of fi,..., and the overline denotes 
the closure in Lf. 

Now, we formulate our design problem. 















































Fig. 7. Discrete-time eiTor system 

Problem 1: Find the digital controller (canceler) K{z) 
that stabilizes the feedback system in Fig. |6] and minimizes 
the following cost function: 

J= sup llzll. 

This is a sampled-data control problem. In the next 
subsection, we reduce this problem as a problem of discrete¬ 
time ^^-induced norm minimization. 

B. Transform into Discrete-Time Control Problem 
As discussed above, the received signal w is in the 
subspace W. Since the delay mT is an integer multiple of 
T, we can represent the delayed input w{t — mT) as 

OO 

w{t - mT) = '^Wd[n](j)n{t), (2) 

n—0 

where Wdin] € is the n-th coefficient. Also, since the 
transmitted signal u is the output of TL^, u is also in W. 
Thus, we have 

CO 

u(t) = ^ Ud[n](j)n{t), 
n—0 

where Ud\n] € is the n-th coefficient. Define 

Zd[n] = Wd[n] - Ud[n\, n = 0,1, 2,... (3) 

Then we have 

OO 

z{t) = w{t - mT) - u{t) = ^ Zd[n]4>n{t)- 

n—0 

Now we have the following theorem: 

Theorem 1: The sampled-data error system shown in 
Fig. |6] is equivalent to the error system shown in Fig. |7] in 
the sense that 

J = sup ||z|| = sup llzdllr^. 

W&W,\\w\\ = l ,\\wd\\i2=l 

Proof: Since W C is a closed subspace, W is 
a Flilbert space. Since spanj^O; </'!; • ■ •} is dense in W, 
{(pOjfi,...} is a complete orthonormal system in W. Thus, 
for any w G W, Perseval’s theorem [8] gives 

\\M\ = \\wd\\e^, 

where Wd[n] = Also, since z G W for any w G W, 

we have 

ll^ll = \\zd\\e2, 





Fig. 8. The band-limiting and pulse-shaping system 'Hh 


where Zd[n] = {z^fn)- Then, w G kF is a necessary and 
sufficient condition for Wd G and hence we have 

sup || 2 ;|| = sup \\zd\\p- 

wClW,\\w\\ = 1 Wd^P ,\\wd\\f2=l 


From this theorem. Problem 1 is transformed into the 
following discrete-time problem: 

Problem!: Find the digital controller (canceler) K{z) 
that stabilizes the feedback system in Fig. |7] and minimizes 

Jd= sup \\zd\\p- 

Wd^t?,\\wd\\i2=^ 


IV. Implementation-Aware Design 


In this section, we will give a design formula for the 
optimal digital canceler K{z). At the same time, we take 
account of implementation of the generalized sampler 5^ 
and the generalized hold 

In practice, the ideal uniform sampler with an anti-aliasing 
filter is widely used in digital systems in place of the 
generalized sampler S^. The ideal (2-dimensional) uniform 
sampler Sh with sampling time h > 0 is defined by 




yi{nh) 1 

) n=0.1,2,...' 


We assume that the sampling time h is determined by 
T = Nil for some natural number N, that is, h = T/N. 
We denote the anti-aliasing filter by F{s), which is assumed 
to be a 2-input/2-output, proper, stable transfer function 
matrix. That is, the generalized sampler is implemented as 
ShF{s). On the other hand, the generalized hold is also 
implemented by using an upsampler f N with upsampling 
ratio N, an FIR (finite impulse response) digital filter P{z), 
and a zero-order hold TTh with sampling time h, as shown 
in Fig. 0 Here we define the upsampler t W by 


t N : {a:[fc]}^o {a:[0], 0,..., 0, 0,...}, 


AT-l 


and the zero-order hold Hh by 


Uh : 


Udi [n\ 
Ud 2 [n] 




n=0,l,2,, 


Ui{f) 

U2(t) 


t^[0,oo) 


We denote the system shown in Fig. |8]by Hh, that is, 

nh = nhPiz)it N). 


This structure is widely used for pulse-shaping [9], and the 
filter P{z) is called a band-limiting filter. The filter P{z) is 


































Fig. 9. Basis function 4> with finite support 



Fig. 10. Discrete-time eiTor system with 'Hh 


usually designed such that its impulse response is equal to 
the sampled values of the baseband pulse cj). More precisely, 
we set 

Piz) = ^(l>{nh)z-\ 

n 

We assume that the support of (^{t) is finite and the support 
length is an integer multiple of h, that is, there exists a 
natural number Np such that the support length of (j>(t) is 
Nph. See an illustration of (j){t) shown in Fig. |9] By this 
assumption, the filter P{z} becomes an FIR filter. Fig. [TO] 
shows the system in which the generalized sampler and 
the generalized hold TL^ are replaced by ShF{s) and TLh, 
respectively. 

Then, let us consider the generalized sampler for 
the error 2 (see Fig. [T]). This sampler is virtual (i.e., not 
implemented), and we can adopt more precise approximation 
than ShF{s). We here assume that the pulse (j){t) with 
support length Nph has its peak at f = Nph/2, the center of 
the support set (see Fig. |9]l. Then, we can well approximate 
S^z by z{Nph/2). That is, we approximate by Sh = 
We show the corresponding block diagram 
in Fig. [To] where T ^ is a downsampler defined as 

4 . N : {x[A:]} I—>■ {a;[0], a;[iV], x[2iV],...}. 


system is defined by 


lift 


■ A 

B ' 

C 

D 




= L^ 


■ A 

B ' 

C 

D 


T -1 

I 


an 

AN-ib 

A^-'^B . 

. B ' 

C 

D 

0 

. 0 

CA 

CB 

D 





. 0 

CA^-^ 

CA^-^B 

CA^-^B . 

■ D . 


where 


LN = {iN)[ 1 2 


= [ 1 




2 

r-N+l 


N-1 iT 


2 -" 2 ---^ (tiV). 


Then we have the generalized plant representation shown in 
Fig. [TT] where we define 


S = 


Til Ti2 
T21 T22 

-1 


K{z) ^ TnK{z)1.-^ , 

Til ^ 

Ti 2 = -SNhe^^’^/^^nhP{z){^N){iN)-Lj,\ 

T 21 = hNShF{.s)nhP{z){^N), 

T22 = LNShFis)Qis)nhPiz)itN)iiN)L],\ 

Note that if 

—- < mN, 

2 

then the system E becomes a linear time-invariant discrete¬ 
time system. Thus, by using the generalized form, one can 
obtain the optimal controller K{z) to minimize the F[°° 
norm of 

via a standard discrete-time H°° optimal control method 

[ 10 ]. 


This system has two sampling rates and becomes a multi-rate 
system, which can be equivalently transformed into a single¬ 
rate system by employing the technique called as discrete¬ 
time lifting [10]. The discrete-time lifting of a discrete-time 


V. Simulation 

We here show a simulation result of self-interference 
cancelation to illustrate the effectiveness of the proposed 
method. 



























































Fig. 12. Impulse response of the band-limiting filter P{z) 


We chose the following simulation parameters: the sam¬ 
pling period h to be normalized to 1, the normalized carrier 
frequency / is 10000 [Hz], the symbol period T = Nh = 2 
(the upsampling ratio N = 2), and the process delay m = 5. 
We assume the anti-alias analog filter is given by 


and the transmission gain a = 2500. We also assume 
the multipath number M is 2, and its characteristic is 
0.2e“^°® -f The baseband modulation is taken 

as the orthogonal frequency division multiplexing with the 
binary phase shift keying [11], where we chose the block 
length as 64 and the guard length as 16. 

We employ the band-limiting filter P{z) as 

P{z) = 0.0165 - 0.0533Z-1 - 0.0177^-2 -p 0.08212-3 
-f 0.01862“^ - 0.14552-3 - 0.01922-3 -f 0.44992-'^ 

-f 0.72872-® -f 0.44992-3 - 0.01922-3° - 0.14552-33 
-f 0.01862-32 -p 0.08212-33 - 0.01772-34 - 0.05332-33 
-P 0.01652-3°. 

This filter is a SRRC band-limiting filter [12] with roll-off 
factor P = 0.1. Fig. [T 2 I shows the impulse response of P{z). 

Fig. [ 13 ] shows the absolute value of the error signal z{t) 
with a log-scale vertical axis. The dash-dot line represents 
| 2 (f)| by the canceller designed by a conventional method 
proposed in [5], and the solid line represents | 2 (f)| by the 
proposed canceller. This figure shows that the error is much 
reduced with the proposed method than with the conventional 
method. From this result, our method that takes account of 
the signal baseband subspace W is much more effective than 
the conventional method. 

VI. Conclusion 

In this article, we have proposed a new method for self¬ 
interference cancelation taking account of the baseband sig¬ 
nal subspace. The problem is first formulated as a sampled- 
data H°° control problem with a generalized sampler and 
a generalized hold, which can be reduced to a discrete-time 



Fig. 13. Absolute value of en'or signal | 2 :(f)|: with Sampled-Data Method 
(dash-dot line) and Proposed Method (sold-line) 

£2-induced norm minimization problem. Taking account of 
the implementation of the generalized sampler and hold, 
we adopt the filter-sampler structure for the generalized 
sampler, and the uspampler-filter-hold structure for the gen¬ 
eralized hold. Under these implementation constraints, we 
reformulate the problem as a standard discrete-time H°° 
control problem by using the discrete-time lifting technique. 
A simulation result is shown to illustrate the effectiveness of 
the proposed method. 
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